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Abstract 

I give a simple proof that the correlation functions of many-fermion sys- 
tems have a convergent functional Grassmann integral representation, and 
use this representation to show that the cumulants of fermionic quantum 
statistical mechanics satisfy i!^-clustering estimates. 

1 Introduction 

Formal coherent -state functional integrals of the type 

Z = Tre-''^^-^^) (1) 

= / n d^(T,a:)d0(T,a;)e-^''^^(^(^'^)(*+^)^(^'-)-^(^(^)'*(^))) 

Te[o,/3) 

have become ubiquitous in the physics literature. These formulas, and their 
generalizations to correlation functions, have become powerful heuristic and 
calculational tools in theoretical physics. However, even when space A is re- 
placed by a finite set, their mathematical meaning is not obvious because the 
product measure over a continuum of r does not exist. The natural way to 
obtain a mathematically well-defined identity is to use the Lie-Trotter product 
formula to introduce a discrete Euclidian time r and recover a rigorous vari- 
ant of the above formal equation in the limit where the time discretization Ar 
vanishes. This involves the solution of a mild, but not completely trivial, ultra- 
violet problem, which arises because only a first-order derivative with respect 
to r appears in the exponent, so that the covariance of the Gaussian integral on 
the right hand side is not absolutely summable as a function of the frequency ui 
dual to T. This problem has been solved both for fermions and bosons. The so- 
lution for bosons is very recent [5]; for fermions, a solution that does not require 
rather obscure multiscale arguments was first given in [4] at zero temperature 
(/3 oo) and then in [3] for all /3 > (the basic problems and the background 
are discussed in [3]). 



The main technical difficulty in the study of condensed matter systems is, 
of course, not this ultraviolet problem, but the infrared problem associated to 
perturbations of Hamiltonians that do not have a gap above the ground state, 
entailing a slow, non-summable, decay of the above mentioned Gaussian covari- 
ance at long times and distances, and much effort has gone into understanding 
its physical and mathematical consequences. On the other hand, proving simple 
statements aboiit the reduced density matrices of equilibrium states requires a 
solution of the UV problem as well. This may be the reason why no simple 
functional integral proof of decay properties of truncated correlations (a.k.a 
cumulants) seems to have been published. 

In this paper I prove the convergence of the above functional integral repre- 
sentation for fermions on a finite lattice. The bounds for the generating function 
for the reduced density matrices will be uniform in the lattice size, so that the 
thermodynamic limit can be taken. 

This paper is partly motivated by the work of J. Lukkarinen and H. Spohn, 
who control the many-body time evolution on the short kinetic timcscalc [6] 
under the hypothesis of an £^ clustering property of the truncated expectation 
values. The result of [6] only requires estimates at some fixed inverse tempera- 
ture f3, but not sharp bounds on the dependence on /?, so that, using the results 
of [3], multiscale techniques can be avoided altogether. Inspection of the hy- 
pothesis of [6] reveals that the clustering property needed there is a special case 
of finiteness of the natural norms used in the analysis of fermionic systems in 
[4] and [2, 3]. Thus the results given here are easy applications of the results 
of [2] and [3], and the main purpose of this paper is therefore to make them 
more widely accessible. To this end, I also include in the Appendix a largely 
streamlined derivation of the functional integral representation along the lines 
of [7], Appendix B, and a new proof of its convergence using the bounds of [3]. 

2 Setup and main results 

Let X be a finite set. In typical applications, X = Ax {1, . . . , n}, where A is a 

spatial lattice and the second index labels all internal degrees of freedom, e.g. 
spin. This choice is not necessary; one could also think of X as indexing any 
orthonormal basis of a finite-dimensional Hilbert space H, such as the energy 
eigenbasis of a system in a trapping potential (with an ultraviolet cutoff). Let 
J-" = /\7i he the antisymmetric Fock space over H; it is also finite-dimensional. 
Let a+ and a~ denote the standard fermionic creation and annihilation operators 
satisfying canonical anticommutation relations. 

Let p be a density operator of the form p = Z~-^e~^'^^~'^^'^ where H = 
Ho + V and N is the number operator (we also call Ho — fiN = Kq). On 
the finite lattice, all these operators are bounded, and the partition function 
Z = Zx{P,n) = Tr^e-'^(^-''^) is well-defined. 



2 



(provided that Zx 7^ 0, which will be the case in our applications, where H = 
H*). 



2.1 Normal— ordered monomials 

The reduced density matrices in the grand canonical state associated to H, (3 
and /X are 



7m,n(a;i,...,a;„j;yi,...,y„) = ( II 



(3) 



\i=l j=n 



0,H,X 



If the Hamiltonian is f7(l)-charge-invariant, 7m,„ 7^ only if m = n, in which 
case I denote it by 7^. 

The reduced density matrices of the state can be obtained by taking deriva- 
tives of the following generating functional. Let {c^)xex and {c^)xex be Grass- 
mann source fields satisfying for all x,y G X and all u, w G {— 1, 1} 



for all x,y e X, and define 

Z{c-,c+) = Tr 



g-/3(//-MW) g(c+,a+)xe(c-,a-)x 



(4) 
(5) 



where {c'^,a'^)x = J^da; c+a+ (here we use a continuum notation J^dx 
^"^xex where e can be a suitable power of a lattice spacing). Then 



'Im {'^1 5 * • • i '^m ! 2/1 i • • * 5 2/m) 



Z(e~,c+) 



(6) 



<5 ^ g-1 d 

Grassmann- valued sources. 



with = e ^ • As will become clear in the following, it is natural to take 



2.2 General monomials 

In applications, one often wants to get information about the expectation values 
of monomials where annihilation operators can be to the left of creation oper- 
ators. This easily included in the present setup as follows. Although the two 
exponentials under the trace in (5) do not commute, they are easily combined 
into a single one, as follows. Because the c+ and c~ are Grassmann variables, 

[c+a+,c-a-] = c+c-{a+,a-} = c+c-Sx{x,y) (7) 

which commutes both with {c^,a'^)x and {c~,a~)x- The Baker-Campbell- 
Hausdorff formula then implies that 

g(c+,a+)x g(c-,a-)x = g(c+,a+)x + (c- ,a-)x + | (c+,c-)x 

= e^'^"'"")-^ e^'^""'""')^ e^'^^''^")^ (8) 
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Thus commuting two such exponentials of source terms produces an additional 

factor e"*^^ . 

2m 

Consider now the expectation value of a monomial Y[ ^^'^'(•^i); where a = 

i=l 

(ui, . . .02m) G 1}^™ is an arbitrary sequence of creation/annihilation in- 
dices. For i e {1, . . . , 2m} with cTj = + let Li{(j) = {j <i -. aj = —1} be the set 
of labels of annihilation operators to the left of i; if cTj = —1, set Li = 0. Then 

by (8) 

2m 

It is an elementary exercise to show that for quasifree states, the prefactor 
gives the familiar 1 — //3 structure for the expectation values of unordered 
monomials in quasifree states (used in the context of deriving the quantum 
Boltzmann equation in [8] and [6]). Here fjs is the Fermi function (see below). 



-T T 
= e ^« ^ 



(c. ,cf)x 



2m 



(9) 



2.3 Truncated expectations 

The truncated expectations (also called cumulants, or connected parts, of the 
reduced density matrices), are the derivatives of the logarithm of Z(c~, c+): let 

F(c-,c+) = logZ(c-,c+) (10) 

then 

jlixu...,Xm;Xm+l,...,X2m) = Sc+'Zci^ Sc^^Jl^^^, ^(^"' L+=o c-=0- 

(11) 

By straightforward expansion of the exponential in Z = and evaluation of 
the Grassmann derivatives at zero sources c~ and c"*", one can verify that the 
7^ indeed coincide with the truncated expectations, as defined in [1]. 

Equation (9) makes evident that for m > 2, the thus defined truncated expec- 
tations 7^ are indeed totally antisymmetric under permutations of xi, . . . , X2m 
(as stated in Theorem A.l of [6]): when taking the expectation value of (9) 
and then the logarithm, the cr-dependent prefactor gives a quadratic contribu- 
tion which drops out in the derivatives with respect to the sources c~ and c"*" 
for m > 1. Antisymmetry follows because the product of exponentials is in- 
variant under permutations, and because Grassmann derivatives anticommute. 
Moreover, (9) implies that the truncated expectation values of order m > 2 of 

/2m \ 

the sequence of monomials i JJ (^"^ i^i) ) are in absolute value independent 

\i=i J men 

of the ordering, therefore bounds on the truncated reduced m-particle density 
matrices imply the same bounds for all truncated expectations (for m > 2). 
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2.4 ^^-clustering 

Let Hq be selfadjoint so that Kq = Hq — = [a^ ,£a~)x with a hermitian 
matrix £. Let 

C(r, E) = -l,>o//3(-i^)e-"^ + lr<ofp{E)e-^^ (12) 
where f 13(E) = (1 + e''^)"^ is the Fermi function, and let 

Cr,r'{£) =<C{t-t',£). (13) 

be the standard EucUdian-time-ordered free fermion two-point ftmction. C is 

an X X X matrix with X = [0,/3) x X. Let (5 > be a determinant boimd for 
C, as defined in Definition 1.2 of [3], namely 6 is such that for aU n G N and all 

Xi,...,Xn,yi,---,yn^X 



sup 

pi,...,p„,gi,...,g„eB 



det((p,,«,)C,,,J^<..<„ <52" (14) 

(here B = e C" : W^W^ < 1}), and let 

a=sup t dT [ dx\C{T,£)\^^y (15) 
vex Jo Jx 

be the decay constant of C (see Section 4.2 of [3]). Both of these constants are 
finite for any ,6 > and finite index set X; S is uniform in 6 and in X if fp{±£) 
is trace class uniformly in X (see below) , and a is uniform in /? if C decays fast 
enough (for sufficient conditions for this in terms of £, see [3] and below). 

The interaction V can be a generic charge-invariant multibody interaction, 
assumed to be given in normal ordered form 

/2m 
Y[ dxi v„{xi,..., X2m) a+^... at^a-^^^ . . . a" ^ (16) 

For fm : X^^ ^ C, let 

/2m 
Yl'^^i\fm{xi,...,X2m)\- (17) 
i = l 

For h > and a sequence / = (/m)meN let 

ll/L=El/™ll,oo/^"" (18) 

m>l 

Theorem 1 Let uj = 2aS^^ , and assume that the interaction V has the prop- 
erties V = V* and uj\\V\\.,^g < 1/2. Then for all m> 2, 7^ obeys -clustering, 

%. 6. 

|7^|,^<2(m!)2a2- 5-2-11^113,, (19) 
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and the same bound holds for truncated expectation values of unordered mono- 
mials of degree 2m > 4. For m = 1 (and denoting at V = by jf^^) 

bi-<o\,,oo^^'''^~'ns- (20) 

Proof: Eqs. (19) and (20) are proven in Section 3. That the truncated 
expectation values of unordered monomials of degree 2m > 4 satisfy the same 
bounds as 7^ was shown in Section 2.3. ■ 



Remarks. 

1. For a two-body interaction V = J^^da; J^dy : nx v{x,y)ny :, = 

^ is small for small (and summable) v. Thus, when measuring the 
strength of the interaction, one can regard the coefficient of /i^ as the 
coupling constant, A = It'll For general V, I take the convention that, 
up to the lowest power of h that appears in it, \\V\\f^ is the "coupling 
constant" . With this convention, 

• Equation (19) imphcs (6.19) of [6] because (m!)^ < (2m)!. 

• Equation (20) implies (6.20) of [6]. 

The factorial in (19) reflects the property that the radius of analyticity of 
F{c~,c'^) as a function of (c~,c+) is finite, while the partition function 
Z{c~,c^) itself is an exponentially bounded entire function of (c~,c+). 
The advantage of using truncated functions is, of course, that the radius 
of convergence for 7^ is uniform in \X\ and in N, and that the expan- 
sion coefficients converge in the limits |X| — > 00 and N ^ 00, while the 
expansion coefficients in Z diverge for jXj ^ 00. 

2. Similar bounds hold for states that arc not charge invariant; the restriction 
was imposed here mainly for notational simplicity. 

3. Estimates for the determinant bound and the decay constant are given in 
[3] for the translation-invariant case where A is a discrete torus of side- 
length L, so that the propagator C has a standard Fourier representation. 
In particular, in this case, S is uniform in L. The decay constant a will 
in general depend on (3, but there are a few cases where it does not, in 
particular if there is an energy gap, i.e. the distance of to the spectrum 
of £ is bounded below by a fixed positive number. 

4. More generally, let {ifix)xex be an orthonormal basis of H, and let 

C(r,x),(r',x') = {Vx\C{t - T',£)iPx'). (21) 

The spectral representation £ = ez|ei)(e;| then gives a natural Gram 
representation (with e"'*"'^ replaced by (e; | ipx) and integration over p 
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replaced by summation over I) that replaces the Fourier representation in 
Lemma 4.1 of [3]. Proceeding as in the proof of this Lemma, one obtains 

5<2max^/0(aeO (22) 

which is bounded uniformly in X if /^(±f) is trace class uniformly in 
X. This effectively requires an ultraviolet cutoff because ff){—E) 1 as 
E ^ oo. The decay constant depends on the properties of the ei and, 
again, on the presence or absence of spectrum of £ near E — 0. 

5. More detailed bounds can be obtained by the methods described in Section 
3: by using weighted norms, one can obtain precise decay estimates; by 
using more details of tree expansions [2] one can also obtain pointwise 
estimates. All these estimates depend on the decay constant and the 
determinant bound of the covariance. In multiscale approaches to the 
low-temperature infrared problem, the bounds given in Section 3 recover 
the full power counting, as already noted in [3]. 



3 Functional integral representation 

To analyse connected functions, it is convenient to have a graded algebra in 
which all even elements commute. This is never the case in the CAR algebra 
of fermionic operators on Fock space. This is the main motivation for using a 
Grassmann integral representation of the grand canonical traces. 
Z{c^,c^) has the following Grassmann integral representation. 

Theorem 2 Let Hq be self adjoint so that Kq = Hq — nN = [a^ ,£a~)x with a 
hermitian matrix £. Let V be any linear map from T to T . 

For m e {1, . . . , iV}, let = and be the operator on (g) C'^l 
given by 

R'£l=Cr„,r^{£) (23) 

withC given by (13). Then 

Z{c-,c+)=Zoe^^''M^'^^^^^ lim (/x^w * e"^) (^W, r?W) (24) 

AT— >oo 

where for G {1, . . . , A^} 

4'^^=Cj,..c^ ,r=C.,,.,c+, (25) 

(the T denotes the transpose as an operator on C'^^l ), Zq = det(l + e~^^) = Tr 
Q-0{Ho-ptN) partition function for V = Q, and 
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Here M{V) is the normal ordered form of V, and * denotes the standard Grass- 
mann Gaussian convolution integral (see [7]) 

(Mc * F){1 = y dMc(V^, + + (27) 

where the integration is over the Grassmann variables (V'fe.x) V'fe,x)(fe,x)6{l,...,Af}xx• 



Theorem 2 is proven in the appendix. 

In the Theorem, we recognize e~^*"'(^'"''^) = /Xjj(iv) *e~^ as Wilson's effective 
action, which generates the connected amputated Green function of the time- 
discretized theory. In [2] and [3] , this function was shown to converge uniformly 
in X, provided that S and a are uniform in X. The variables r/^^^ and r]^^^ on 
which W depends contain a C because Z generates non-amputated functions, 
where a propagator is associated to each external leg. The sources c appear in 
the time slices ri and tn = as expected since the monomial "sits" at the 
fixed time r = 0. 

I briefly recall some of the properties of W proven in [2, 3] . By Theorem 4.5 of 
[3], the limit W(C, V) = limjv^oo W^^\R'^^\V) exists in \W for small enough 
h. Let h = 6 in that theorem and denote the part of W that is homogeneous of 
degree p in V by W{C,V;p), then 



P 1 

W(C,V)-^->V(C,V;p) 



1 P- 



l-'-IIVIIs.x'""^^'^ 



< ; "..i^r mss.x (28) 



holds for all P > 0, provided t^HVUg^x < i^*^^ P = the empty sum in (28) is 
zero; this is the only case needed to prove Theorem 1). 

The norm |Hl/jx similar to the norm defined in (18), but with ^ 
in (17) replaced by a l,oo-norm Hi^^x which integrals and suprema run 
over X = [0, /3) X X. Because the interaction V gives rise to a V that is local in 
the timere [0,/3), llV||,_x = 11^1^. 

In summary, for small enough w||y||3^, W is analytic in V and in the fields. 
In particular, W is continuous in the fields, and 

lim W(^) (i?(^) , V) (77(^) , 7?(^) ) = W(C, V) (r), r?) (29) 

JV— >(x> 

with f]{T) = Cj and 77(r) = C^fiC^. Let Wj^'' be the homogeneous part 
of >V(") of de gree 2m in rj, rj. It generates the connected amputated 2TO-point 
functions W^^'' = j^^2m ■ Then W^^^ is given by a tree expansion 

[2, 3], i.e. an explicit, absolutely and uniformly in N convergent expansion in 
which the limit iV ^ oo can be taken termwise. Denote the limit by Wbm- 
Then, ifw||y||35<l/2, 

\W2mU^,^ < 2<5-^- ||y||3, (30) 
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for all m > 1, and 

Z{c-,c+) = Zo e(^"'/''(^)'="')^ e-^^^'") (31) 

so F = log Z becomes 

F(c-,c+) = logZo + {c-,.fp{£)c+)x-W{fj,v) (32) 
= logZo + ic-,[U{S)-IC0]c+)x -Y,W2miC0,-c-,C.,oc+) 

m>2 

with /C/3 = J dr J dr' Cp,rW2{T,T')Cr'.o- In (32), the full two-point function is 
explicit in the quadratic part of the generating function. The analogue of the 
free fermion formula, obtained by discarding the yV2m for all m > 2, changes 
correspondingly. Under the present assumptions, W is small when ||V||^ is small, 
and hence it is clear that all higher truncated functions are small and ICp is small 
in I'l]^ ^^"^ two-point function (m = 1) is close to the free one. 

Theorem 1 now follows immediately from the above properties of W: Be- 
cause the arguments of W in (32) are fj = C/}^. c~ and rj = C.^C^, 

/. 2rn 
Y[ dTfcdyfe {m\fW2m ((Tl,yi), • • • , {T2m,y2m)) 
k=l 

rn 

n(C(/5'-^))x.,..(C(r.+™,0))^^^^,,^^^ (33) 
fe=l 

for all TO > 2, hence 

< |W^2„^|l,oo,X |C|i:^,x = \W2m\,,^,^ a'"^ (34) 

by (15). By (30), 

< 2(to!)2 ||F||3,a2™ (35) 

which proves (19). For to = 1, recall that 7^g = so 7^ — 7^3 = /C/3, as 

defined after (32). Thus, again by (30) and (15) 

- ^mIi,oo = \'^<oo < «'l^2ll,oo,X < ^C.'S-^V\\,„ (36) 

which proves (20). 

A Traces and Grassmann Integrals 

Lemma 3 Let be a Hilbert space and A and B he hounded linear operators 
on T . 

e^+^= lim 

JV— »oo 



e~ (1 + ^) 



(37) 
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Proof: This is an easy variant of the standard proof of the Lie product 
formula, as given, e.g. in [7]. Let C = e(^+^)/^ and D = e^/^{l + B/N). Then 

max{||(7||,||£)||} <e^^ (38) 
and by following the proof given in [7] , one obtains the estimate 

- < ^{\\A\\ + \\B\\f ell^ll+ll^ll. (39) 



By assumption, the one-particle Hilbert space is finite-dimensional, i.e. has an 
ONB indexed by a finite set X. Thus the Fock space is finite-dimensional, 
too, hence all linear operators on T are bounded and the trace is a continuous 
linear functional. Therefore by Lemma 3, 

Z{c-,c+)= lim Zn{c-,c+) (40) 

where 

Zjv(c-,c+)= Tr Ljv e('=^'"+)e('="'")] (41) 

and pN = (CnDn)^ with Cn = e-#^" and Dat = 1 - . If V is normal 
ordered, the same is true for Dj^ (but not for This is the main reason 

that makes Lemma 3 convenient in the transition to the Grassmann integral 
done in the following. 

By a purely algebraic transformation, Zjv(c~,c+) is represented as a finite- 
dimensional Grassmann integral. Some basic facts needed for this are in Ap- 
pendix B of [7]. Here I collect only the most important definitions and identi- 
ties. The Grassmann symbol G{A) of an operator A with normal ordered form 
J\f{A){a'^ ,a~) is defined as 

a(^)(^+,^-) =e('^+'*-)- AA(A)(^+,V-). (42) 

The product of two operators has Grassmann symbol 

a(AB)(V^,V) = J I'x(^', V) Q{A){il',^') e-^^''^')^ G{B){i,',i,) (43) 

The trace is represented by 

T^^A = J g{A){-i;+,i,_) (44) 

By these identities, Z can be written as the convolution 

Zn{c-,c+) = J vx{0,e) g{pN){e,e) e^s+c- ,e+c^) ^ ^^^^ 
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By the rule for products of Grassmann symbols, (43), 

r N _ 
•' fe=i 

N N-l 

\{g{CNDN){^k,i^k) [] e-('^'=+i''^'=)- (46) 

The decomposition H = Hq + V ]s chosen such that Hq is quadratic, so that 

K{) = {a'^ ,£a^)x with a licrmitian matrix £. For this special case, the Grass- 
mann symbol is easily obtained by Grassmann Gaussian integration, with the 
result 

g((7;v£'iv)(4, V'fe) = Q{DN){uN^iJk, A) (47) 

where Un = e~^^ (and T denotes the transpose). The interaction V is assumed 
to be in normal ordered form, so its Grassmann symbol is straightforwardly ob- 
tained by replacing a"*" by ip and by ip. Inserting this and changing variables 
from -iph to Uisr'^'ipk gives 

r ^ - 

ZM{c-,a) = dete-M []l?x(i^fc,^fc)e("~"*+'="'*"+'=^)^ 

k=i 

^ / (J _ \ w-i 



fc=i ^ ' fc=i 



N 



e(="'^^)- dete-^ / flVx{^k,^k) e^'^''^' 
•' fc=i 

n (l - ^V{{pu, i^k)) e(-"'^«)-+*'«"-+)- (48) 
I — 1 ^ / 



with 



fe=i 



JV N-l 



(V',(3^^V) = '^{4'k,'4'k)x - ^ (V'fe-M, wjvV'fc)^ + {■4>i,UNtpN)x (49) 



fe=l k=l 



That V^i couples to V'jv with the opposite sign means that there is an antiperiodic 
boundary condition for the fields ■i^ and ip. Q^^^ acts as a matrix on the "time" 
indices k and as an operator on the x £ X: 

Qm}n = ^rn,n " UN^m-l^n + WAT^m.l^AT.n = Qm'^^ + UNSm,lSN,n (50) 

q(o.^) is lower triangular and its inverse is easily found to be it!^^'^' 

4°'^) =«jv'=-' ]l(fc>/) (51) 
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The additional term in Q is similar to a rank one perturbation. Using the 
standard formulas, one obtains {Q^^^)~^ = R^^"^ as 



(ll(m>n) + 

— {1+un~^)~'^ for to < n 
{l + UN^y^ forTO>n 



= Un'' 

we obtain the formula (23) for R^^\ Thus, in summary 



(52) 



diiR(N){i),il)) e^' 



fe=i 



N 



(53) 



where Z^^^ = dete"^ detQ^^^. 

It remains to reexponentiatc V. To this end, recall from Definition B.ll of [7] 
the £^ seminorm |||-|||q on the Grassmann algebra, namely = '^rn>i l^mli<?™ 
with \Am\i = JYl^i<ixi \Am{xi, . . . ,Xm)\- It satisfies the product inequality 
< III A III g III B III g, and if the covariance C has determinant bound 6, then 



(see Theorem B.14 in [7]). Now let 



< IIIFII 



AT 



TV 



N 



A(vl/) = []e-^^(*^)-n 1 
fe=i fe=i ^ 

By the discrete product rule 

N N N 

fc=l k=l 1=1 k<l k>l 

and by (54) and the triangle and product inequalities, 



(54) 



(55) 



(56) 



/ 



dMc(*)A(*) 



< 



IIAIIIi 

(JV-i) 



< iVe^/'ll^ll^|||e-^^-l + ^y 



(N-l) , 



, nv\\ 

7V2 



e" ' 



N 



1 e^mu ^ 



(57) 



as ^ cx). Note that if X is a lattice of volume with a translation group, 
III G III g is typically of order Q for translation invariant elements G of the Grass- 
mann algebra, but this does not matter here because O remains fixed as — > oo. 
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Thus 

Z{c-,c+)= lim Zn{c~,c+) (58) 

N—^oo 

where 

Zn{c-,c+) = e^^"''^^)^ det(l + e-^) y d/^^w (^i, V) e('^"''^")^+('^^''^^)^ 
e-#ES=iV(v;fe,Vfc) (59) 

with given in (23). Eq. (24) now follows by a standard completion of the 
square in the integrand. 
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